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Abstract 

ouallti« n ri?J2HJl e 2 p ??" lo ! p *Pf c1 ff cations, for example, aircraft handling 

des 1 2 n *«****«. particularly modern control approaches, 
SiT,;: a 11 a * bl V° t 5 e systeffl signer for synthesizing even the most complex 
4 contr ?] systems. Unfortunately, however, weaknesses exist in the 

If ?n 9 i q :: itfa ! a , rea$ ‘ particularly for "non-conventlonal" aircraft such 
as V/STOL and control configured vehicles (CCV's). In this paper, an auo- 

Qualities sMJf«J%T eth0d u ! abl ! in f he absence of quantitative handling 
SlU*!! JF l J 1f1 f at1ons » and y«t explicitly Including design objectives 
based on pilot-rating concepts, will be presented. The algorithm Involves 

sistS ?SAS? P !a r ilf Sol ): in ?, for the stab ^ty augmentation 

IS 1 * 9 1 ?* Pilot equalization and pilot rating prediction via optt- 

thl 5 <m ^t an *°us solution Is required In this case since 

thJ (jalos, etc.) depends upon the augmented plant dynamics* and 

the augmentation Is obviously not a priori known. Another special feature 

m*5«) S to' , de^9n ,, {hJ t SAS 0 ‘ >J * CtW * f “ nCt, ° n (f ™” Wh,ch the p,,0 ‘ ITOd *' 


Introduction 

.. 6 | ven a i equat ® design specifications, or aircraft handling qualities 
i 9 anc ^ a va ^ system model, design techniques, particularly modern 
control approaches, are available to the system designer for synthesizing 

nI«J he 7°i t < C °?u 1e 5 f ! ight contro1 systems. Unfortunately, however, weak- 
J^h e L*v/sTfti 1n t 5 e de ! 19 ? s P e *j fi cat1on area for non-conventlonal aircraft 
c h liVt* iEf «° ntro1 con ”« ur ? d vehicles (CCV's). The assertion here 
Li h JL d n?i^<! he * "°"- co "y«nt1onal," multi-variable nature of the vehicle 
tas J 1n f h * $ ase of ^/STOL ) , and due to the anticipated 
?^ y ? f the ^ ste ? s 1nvolved » a non-conventlonal" approach to the 
control design problem is worthy of Investigation. 


„<i„* S1n S e < p11ot acc fP tan « Is the ultimate criteria, In the absence of prior 

P !l!°!! w * mu$ J predict pllot rat1n 9* This 1s in contrast to deslSn 
*nH b ? d | attem P t *° a P r lori define "good" dynamics, and then use a 
model -foil owing design technlquellJ, that Is, design the auonentatlon so the 

thU^M^ch 1 ?! Th 1 ! beba lf e 11ke the “ 90od " mode1, 0ne ma ^ or drawbac k to 
this approach Is that one Is never sure that the pilot will agree with the 

fcOV 

?MC &TENTIONAUY BUI< 
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designers choice of "good" dynamics. 

types of^plfot ^models exist 9 * larh K 0rm P^°t wodel Is required and two 

^d2?«55? y; I r?? x?; ude 

Is more compatible with the muTt1-varTabie°«n <> \ Pll 5 t *II 0<,el 1s 1dea1 It 
advanced control design techniaL« V »i!f*] e as P e f^ s of the problem and the 
the pilots equalization n*Wk?« lll“ d L e * ? tf " 9 * Also * the form of 
tant property In this case S automat1ca11 «f determined, a very Impor- 

closedOoojTsystS^^ ha * e be f? used Primarily to study 

control modeMnclode foM^vtsabl 1 1 C J t A?Si K1 as for the opt1ml m 

and the stability of the pilot 1 illi of CTOt^J and STOl aircraft 1 - 53 

K~4.w&ir i. 

SXT th * perfor- 

th. sS*s5s ^:5 s h * 5 f o; n ^ <, :p^t p :i t .;?.5; rt , of ‘ h % p '*'’' C73 - «• 

the pilot model must be assumed before beolnni™ fM?"!’ < Bwy the form of 

sfjJK’ ?.«: zrsrj: 

of 

optimal-control pilSt lodi?! $ wturally accomplished with the 


to predict plfo^opfnlon.^The mw^notabirorthesft^h 9 ? " 0t as fre< * uent ^» 

the VTOl hover test, was the %Tpl r JllJr dmfoSed j! C f n n lj ueS ’r8 P 9l 1ed t0 
this approach, parameters in thl 2L oeveioped by Anderson l a* 9 -*, in 

are chosen such that a pnot-ra^ng metrie fs 5 !?^ f ‘f nc J 1or ’ ® f assumed form 
slsts of a measure of performance ?^ Ll \ ' 7 1 Zed * T C 1s ■•trie con- 

of ptlot workload the aoouJ J'l'.IJ tKi!!! £T! r li *"! * »«»« 


Of pilot workload (o.o„ tSoii!S.l’J , {.l?.k p, *fJ"? ,rror )' ‘"<1 o wasur 
* * of tih ichM J3) ? t f ° oust Introduce). In 

10 ..ally m£X4mtui 2LSS!e“2 '!«»». f«»‘- 


in assessment 

tlonal based on easi lymeasured l mot1 on ZJSSL'** a r.tlni'W 
opinion prediction. However "titles was adequate for pilot 

some additions for better system narformanr^ was ^ ound to require 

Improvements Included mSlffcetloIrln^fe!? ?nfr^ d J K? 1on i* Motab, y» these 
addition of the pilot ■ • 7r«SSt Sr rt b1 5? funct1on ) a " d the 

trol Input.) Hence, again wJ s*l°fhi h ^r? dw P° r t1on of the pilot's con- 
form of the pilot's describing function** ** Created by <m P° s1n 9 an assumed 




- csn^A'ai swsuw?aw»rs?ag£, 
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In analytical display design for helicopters^ use of the nnHmai 

Jl!i, tn f* p lI ot ?°:*' °W ect ^ ve function. Proper selection, based on the 

,*?& ?orh.1£™?l» ntro i r; ,9h ‘; ,n «• p?»- 

.* ^ wtcnnlftjtion of the pilot Qd1n$« Gou^llzatlon anrt n^iAf 

rating prediction simultaneously. As we have mentioned previously 9 this 

dynanlcs^for^hl ch^he^ nl lot ' m ! ystems w1th non -conventional 
oynomics ror which the pilot's describing function may not be known if 

sygr* - - ,th thes « o* 5 ' 9 " 


The Pilot Model 

assumDtlSJ» e thJt e thl n w5ii fe r! n f e 2' th ® °P t1mal Pilot model evolves from the 
InMtfe! 0 th well-motivated pilot selects his control 

,0 h “" ,n ,, " 1t,t,ons > s»cl> th.t the following ob- 


i T 

J p * 1 T / <y'Qy ♦ 0p«Op + o*Gfi p )dt} 


1 w"reUtlS teW be1 " 9 controned by the P ilot is Ascribed by the familiar 


X 


■ V + Vp * 5 


y • Cx 


0) 


where it Is the system state vector 
put vector, and 9 Is the vector of 
variance 


• iip the pilot control vector, y the out- 
zero-mean external disturbances with co- 


ECw(t)w'(t+o)J - H$( 0 ) 


Included as human limitations are observation delay, x, and observation noise 
0 f So the ptlot ® ct «®Hy perceives the noise-contaminated, delayed states!*’ 


y p * C p x(t-x) ♦ v y (t*x) 

of Is^denoted* ^ *"**• 


E[v y (t)v y (t+o)] ■ V y «(o) 


Oeflnlng the augmented state vector, 
or the pilot's control is given as 


X • col[x,u p 3, 


the solution to the problem. 
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5 p * ■ -s''[»;i]y 

where K p Is the positive definite solution to the Rlccatl equation 


TAn J B“1 

1 

K - K 

\k n ! 

P ! p 


V " C P ! 5 


P P 

_° : o 


.... - 


KBG^B'K = 
p o op 


It will be convenient to partition K 



P 


such that 


( 2 ) 


| 

i 


i 


and note that now the equations for the optimal control u p * Is 

u* « -G _1 K n £ - G"V u* + 

P P 3 P 4 P 


or a linear feedback of the best estimate of the state, x, and some control 
dynamics. (These control dynamics have been shown to be equivalent to the 
pilot's neuro-muscular lag.) 


Now, the state estimator consists of a Kalman filter and a least-mean 
square predictor, or 



X(t-t) « A p x(t-T) + sCp'V^CYpft) 

- CpX(t-T)] + BpUj(t-T) 

x(t) - B(t) + A T [x(t-t) - B(t--r)] 


tTo model the pilot's remnant, motor noise Is usually added to the control 
equation. The final pilot's control Is represented by 


i 

u 


* 

P 





K 





where 


i 

! 


ECvJtJVm'U+o)] - V m 6(o) 
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if.* wu V 


•V 


, 




with 0 - AJ + B u p 

+ r is the solution of the Riccati 
,„d the estimation error covariance matrix £ Is 

equatl on 

......... - — ■ *— “ “ 

gaSSSSri^i sasrttw s 

Cooper-Harper system, the relation Is 


Pilot Rating (PR) * 2.53 In (10 Jp) + 0.28 
the dynamic system. 

Au gmentation Synt hesis Method 

In the determination of the ^trlce^An andlp .^However, 

MK augmented pAt. *p and 

B D Is as yet unknown. 

Consider the un-augmented plant dynmnics to be described by 
j = Ax + Bu + w 

where, as before, x Is the system state^vector^and syl^matrlces, and 0 is 
vector. However, A and B are now^t contro i input to the plant will 

SSdTXt \mV C5’“ 4U ^ ntat, °" 1np,,t • ^ or 

U - Up + G SAS 

Further, from th^pllot mode ^^e^euml^fhe pi ?o?f control input 

(e.g., G K p,» u 

In expressible as 

m a _1 (3) 

Up * ‘ 6 ’ K p 3 * ‘ G V p 


*«., r-AH. ‘ WJhrtW'AAWn 





Now, the estimate of the state Z ~ >. 

state plus some estimation error? ’i^J P be expressed in terms of the true 

A 

X * X + e 

control equation * S another disturbance, w p , we can write the pilot's 
U P * ~ 6 " 1 k p 3 * " Up + w 


• 

i 

t CD 
1 


B~ 


1 

X = 

rs ! Aj 

'itj 

_°J 

U SAS + 

-V 


( 4 ) 


where x * col[x,u 1. 

P 


W* n0W " ay Pr ° Ceed “ obJect <ve function f or d ete™,„,„ g 


Si £S 5 ?x*s& S 

Q and R, as the inverse of the mavimiim »ii L?* abe and control welohts 
as perceived by the pilot ) FiIJii!! Um * a owable deviations in the vari«hi«« 
gains, we must also penalize auamllfi+J 0 preclu de infinite augmentation bl6S 
augmentation is cfc^'SJK^*''" C °" tr< " “«W- ThI?l?o?e!°?he 

J SAS * Jp ♦ E (1 Im f / J' As F J sas d{) 

J SAS ■ E Ota T ! 0 <»•« ♦ “p R Sp * UpGCp ♦ 5sA S F=s A s)dt> 

" I"' a " d 6 *" “ *' «• Plot's objective function, 0 . 

This may be written as p 


or 


21 ? 




where 


J SAS = E T f Q ^' P * + u SAS Fu SAS^ dt} 


CpQCp + K^G K p ^, KJ^G K p ^ 

; t :r _ 

Ki g” k p I r+k: g 'k p 

P 4 P 3 ! P 4 P 4 


and Instead of Equation 3 being substituted for u p In the above J$AS* we have 
Invoked a sort of separation principle and substituted the relation 

“p ■ - g '’ k p 3 5 - g \ 5 p 

The justification for this relation being used lies in the fact that we wish 
to synthesize the augmentation based on how the pilot is trying to perform 
the control function rather than on how the pilot is capable of doing so. 

With this objective function and the system dynamics given in Equation 4, 
the problem is now stated in conventional form, except K p and K p are as yet 
undetermined of course. If we assume, for example, full state feedback, 
the solution of this problem is known to be 

“SAS * i 0]K SAS* 


where 


k sas 3 ! k sas 4 


5 sas ■ - f ‘’ b ' k sas/ - F '’ B ' K SASj 5 P 


*sas, ! k sas 


is the solution to the Riccati equation 


- — » — -4— — ?■*” 

-G~K p | -G -1 K p 
P 3 ! P 4 


Sas ■ iAs r::T 


-G" Kp I 
P 3 1 


-G" 1 K, 


p + k sas j* f " 1[b ’ * 0]k sas * *sas 
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on K p (or K D and K D ) P Returnln^t^fh °D? tl0n /° r KsAS obviousl y depends 
ga 1 n P (Equat? 3 n 2 ), Si alsosel Rlc f* ti equation for the pilot 

SAS gains (or KsaO since th* n*w * ?? uation depends in turn on the 

sarrs Ss 0 r d h^eve^ f h ri^ tC ? ? 

cii wrr/ist si? 


A p - A 


BF'Vk, 


SAS, 


and B p = 8(1 - F’Vl^ ) 

COUpl^(«S«l qu t a?o l .s“ P SM S S? S thJ ?nOt'Stn? IC Ff < ff ,Ua ^° n yfe ' dS *» 

the SAS gains, Equation 5 These J,!* h* I?i, 9 j 1n !’ j3 uat 1on 2, and one for 

Kp by Integrating both equations backward 'Note S thatShir?li « r K § AS antl 


is represented in 


A Simpl e Numerical Example 


namtcs°cons1dered*1n l Ref™ll| n9 tSSk WUh * he Controlled dy- 


0(s)/6(s) = K/s‘ 


11.7 


The command signal, e c , is white noise, w, passed through the filter 
M s )Ms) s 3.67/(s 2 + 3s + 2.25) 


a «d E(w)-o, a 2 =1.0 

w 

It we define the state vector as 5 . col(e c ,4 c ,o,e), we have the plant 
X » Ax + B$ + 3.67w 

where 

o 1 0 0“ 

-2.25 -3. 0 0 

A = 

0 0 0 1 

_ 0 0 0 0 


B' ■ [0,0,0,11.7] 


AUGMENTED PLANT 



Figure 1 Piloted Vehicle Schematic 









For this system, error and error rate are perceived by the pilot or 



e - e 


7. 

• 

1 

o 

o 

_J 

«<l 

■O 

u 

c 

_ 9 c ’ 6 _ 

* 

0 

1. 0. 

1 

1 •"* 


The performance Index, chosen consistent with Hess's rating hypothesis Is 

J p = E tllm 1 l [(e - e) 2 + .Q16 2 + gfi 2 ]dt} 
v T-*» 1 0 v 

and g Is chosen to yield a neuromuscular lag, l/x N = 6 ^Kp. = 10., or tjj = 
0.1 seconds. Unagumented, the pilot Rlccatl equations absolved with the 
following noise statistics (human limitations) 

1) Equal attention allocation between error and error rate. 

2) Observation thresholds on error and error rate = 0.5(units of display 

displacement.) 7 

3) Sensor noise, (V y )-M/E(yf ) * -20dB 1=1 ,2 t 

4) Motor noise, (Vu;/E(ug) B -20db, where u c * -G K p x 

5) Observation delay, t * 0.1 seconds J 

The "piloted" system performance Is given In the following table. 


Tabl e 1, Un*augmented System Performance 


(0 C -e) rms rms J p P.R^ 

“TTT" 1.00 U9 


* This pilot rating has been verified by experiment 

Assuming full-state feedback, the augmentation control law is 

S SAS ‘ - K l e c ' K 2®c ' K 3 9 ‘ V ' Vp 
The SAS objective function is 



so the piloted plant, including augmentation will be 
x ■ ApX + BpSp + 3.67w 


where 


0 


0 


0 


A 


P 


-2.25 

0 

-11 .7K 1 


1. 

-3. 

0 

-11.7K 2 


0 

0 

-11.7K 3 


0 

1. 

-11.7K 4 


Bp * [0, 0,0,11. 7(1-K 6 )] 

Solving the pilot and SAS Rlccatl equations simultaneously, and then de- 
1 n™he fol 1 owl ng^ta bl e 6 " 1 perfonnance as before * ields the results given 


Table 2 Augmented System Performance 


f 

(e c - 6) rms 

6p rms 

J P * 

P.R.** 

100.0 

1.10 

0.89 

1.21 

6.6 

10.0 

0.79 

0.61 

0.62 

4.9 

1.0 

▲ _ 1*. a 

0.38 

0.35 

0.15 

1.3 


* not e jsAS S 1S ^ nUI " erical value of the P ilot ' s objective function, J p , 
** Predicted pilot rating based on J p . 


IS%I!E e ! t# !l 0 !! ?*• p - K 4 a " d K 6» for the three cases above are given 

In Table 3, along with the augmented plant eigenvalues. 


Table 3 Augmentation Gains 


f 

i" 

h. 

h 

Hi 

i 

Plant Eigenvalues* 

100.0 

-.009 

CM 

O 

o 

1 

.009 

.003 

.004 

-.017+. 331 j 

10. 

-.078 

-.016 

.084 

.024 

.036 

-.142+. 982 j 

1. 

4 li-x j 

-.513 

-.090 

.542 

.130 

.155 

-. 758+2. 40 j 


* Not Including noise filter eigenvalues of course. 


Summary 

In summary, we have cited the flight-dynamic and control problems of 
non-conventlonal flight vehicles (V/STOL and CCV) due to the compleM v of 
augmentation required and the lack of handling qualities objectives We f 
have presented a methodology Intended to be suitable for this type of pro- 
blem. The method uses an optimal control pilot model, not onlyto predict 
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piloted performance but pilot rating as well. With the optimal -control model 
structure, we were able to formulate the augmentation synthesis problem as 
an optimal control problem with the parameters In plant matrices depending 
on the pilot model, and vice versa. This necessitates simultaneous solution 
of the two (pilot and augmentation) problems. We have Included the form of 
the solution under the assumption of full-state variable feedback and no 
measurement noise, and a simple numerical example. 

.. ™ e first extension to be addressed will be the solution for the case of 

tlt!i« d thJn t ?hI eedbaCk li Th1s . c * se 1s actually closer to pure plant augmen- 
tation than the case addressed here. In our solution, and In the example, we 

nave closed the tracking loop, and pure plant augmentation would only feed 
back plant states. However, the primary purpose of our discussion here was 

to provide the problem structure which would be unchanged regardless of aug- 
mentation approach. 3 

... Further extensions will also include the cases with state estimation, 
with and without measurement noise. Also, the necessity of pre-tuning the 
pilot model will be Investigated. s 
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